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We search for the standard model Higgs boson produced in association with an electroweak vector 
boson in events with no identified charged leptons, large imbalance in transverse momentum, and 
two jets where at least one contains a secondary vertex consistent with the decay of b hadrons. 
Wc use -^1 fb~^ integrated luminosity of pp collisions at ^/s = 1.96 ToV recorded by the CDF II 
experiment at the Tevatron. We find 268 (16) single (double) 6-tagged candidate events, where 
248 ± 43 (14.4 ± 2.7) are expected from standard model background processes. We place 95% 
confidence level upper limits on the Higgs boson production cross section for several Higgs boson 
masses ranging from 110 GeV/c^ to 140 GeV/c^. For a mass of 115 GeV/c^ the observed (expected) 
limit is 20.4 (14.2) times the standard model prediction. 

PACS numbers: 12.60.Jv 
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The Higgs boson is an essential element of the standard 
model (SM) of particles and their interactions explain- 
ing the mass of elementary particles and playing a key 
role in the mechanism of electroweak symmetry breaking. 
Precision measurements of the W boson and top quark 
masses suggest that the mass of the SM Higgs boson 
should be less than 182 GeV/c^ (95% C.L.) jT] and there- 
fore kinematically within reach at the Tevatron. Previous 
searches for the SM Higgs boson set limits on its produc- 
tion cross section at y/s — 1.8 TeV and y/s — 1.96 TeV 
The CDF Collaboration is now performing searches 
for the low mass Higgs boson produced in association 
with an electroweak vector boson V {Z or W) followed 
by the decay H ^ bb in several final states which are 
optimized separately and then combined. In this Letter, 
we describe the first search for the SM Higgs boson with 
a trigger based on an imbalance in energy in the plane 
transverse to the beam sensitive to both ZH and WH as- 
sociated productions. In order to avoid overlap with the 
search for WH in a final state containing a high-energy 
electron or muon candidate, large missing transverse en- 
ergy, and one or two b jets, we veto events with high 
energy isolated leptons. Therefore, we search for VH 
events where H ^ bb and either Z ^ i/D or W ~> f.v, 
where / indicates that the lepton was not identified in 
the detector. The signatures of these decay modes are 
final states with two b jets from the hadronization of the 
b quarks, missing transverse energy, and no isolated lep- 
tons. In this search we use an integrated luminosity of 
0.97±0.06 fb-i collected in Run II of the Tevatron. After 
describing the experiment, we introduce our event selec- 
tion, major backgrounds, and systematic uncertainties, 
and present the results. 

CDF II is a multipurpose detector described in de- 
tail elsewhere [2- Here, we briefiy describe the detector 
components used in this analysis. We use a cylindrical 
coordinate system around the proton beam axis in which 
d and are the polar and azimuthal angles, respectively, 
the pseudorapidity 77 is defined as 77 = — ln(tan(|)). The 
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transverse momentum and energy of a particle are de- 
fined as = psin0 and Et = EsinO, respectively. The 
tracking system consists of a cylindrical open-cell drift 
chamber and silicon microstrip detectors in a 1.4 T mag- 
netic field parallel to the beam axis. The silicon detec- 
tors |1] provide tracking information for jryj < 2 and are 
used to detect collision and decay points. The drift cham- 
ber [5] surrounds the silicon detectors and covers the cen- 
tral rapidity region |?7| < 1. 

The energies of electrons and jets are measured in 
calorimeters which cover the region \r]\ < 3.6 and are 
segmented into towers pointing toward the center of 
the detector. Jets are reconstructed from energy de- 
positions in the calorimeter towers using a jet clus- 
tering cone algorithm [6, with a cone size of radius 
AR = -^Z (A(/))2 -I- (Arj)'^ — 0.4. Corrections are applied 
to account for effects that can cause mismeasurements 
in the jet energy such as non-linear calorimeter response, 
multiple beam interactions, or displacement of the event 
vertex from the nominal z = position [6]. The I^t is 

defined as the magnitude of I^t — — X^i E'^f^ij where fii is 
a unit vector parallel to the transverse component of the 
vector pointing at the i*"^ calorimeter tower, and E'^ is 
the energy therein. Both magnitude and direction of the 
T are recomputed after the jet energies are corrected. 

We use loose electron and muon identification crite- 
ria in order to veto events with at least one high pT 
isolated lepton. Candidate muons are identified with 
Pt > 10 GeV/c isolated tracks which leave energy in the 
calorimeters that is consistent with a minimum ionizing 
particle. A track is called isolated if the total calorimeter 
energy in a 0.4 radius cone around the track is less than 
10% of the track pT- Candidate electrons are clustered 
energy deposits in the electromagnetic calorimeter with 
Et > 10 GeV which have electromagnetic to hadronic 
energy ratio and shower shape compatible with electrons 
and are associated with an isolated track. 

The events used in this search are selected by a three- 
level trigger system. The Level 1 trigger requires I^t > 
25 GeV , where the ifx is determined using calorimeter 
trigger towers with Et > 1 GeV. The Level 2 trigger 
selects events with two jet clusters of Et > 10 GeV, 
where one of the clusters is required to be in the region 
1 77 1 < 1.1. The Level 3 requirement is^T > 35 GeV. After 
selections to remove accelerator-produced and detector- 
related background, as well as cosmic ray events, 4.3 
million events remain. This is about 60% of the initial 
sample. Further selection criteria are defined using cor- 
rected observables in order to search in a region where 
the trigger is fully efficient, but the sensitivity of the 
search is not compromised. We require the events to 
have l^T > 50 GeV and to contain exactly two jets with 
Et > 20 GeV and \ri\ < 2.4 where the azimuthal separa- 
tion of the two jets is A(j>{Ei^T, E2^t) > 1-0 rad. More- 
over one of the two jets must be central with jryj < 0.9, 
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and the jet of highest (second highest) energy in the event 
must have Et,i > 35 GeV {Et,2 > 25 GeV). The trigger 
efficiency for such events is 0.97 ± 0.03 [7]. 

We require at least one jet to be identified as a 6 jet by 
the 5-tagging algorithm SECVTX ^ . This algorithm looks 
for tracks within jets that form a secondary vertex sig- 
nificantly displaced from the primary interaction point. 
The tagging efficiency for b jets of 50 GeV is ~42%, the 
misidentification rate for 50 GeV light quark (it, d, s) 
and gluon jets is ^--^1% |H]. 

In order to avoid potential bias in the search, we test 
our understanding of the background in control regions 
that are defined a priori. The observed events are divided 
into three non-overlapping parts, two control and one 
signal regions, which have different event topology and 
isolated high-p^ lepton multiplicity. The backgrounds in 
the three regions have contributions from multi-jet, tt, 
W plus jets, Z plus jets, and electroweak {WW, WZ, 
and ZZ) processes. Control region 1 (CRl) is defined 
such that it is dominated by the multi-jet background. 
The events in this region do not contain any identified 
leptons, and the separation between the I^t and the sec- 
ondary jet (A</)(i?2,Tj-^T)) is required to be less than 0.4. 
Control region 2 (CR2) is selected by requiring at least 
one lepton and A(j){E2,T,i^T) > 0.4. This region is sen- 
sitive to electroweak and top quark decays. The signal 
region contains events with no identified high-pT^ leptons, 
A(I3{Ei^t,$t) > 0.4, and A(t){E2.T-,$T) > 0-4- 

Agreement between the number of events expected and 
observed due to different processes is presented in Table |l] 
for single and double 6-tagged events in CRl and CR2. 
The signal and the physics backgrounds have been eval- 
uated using PYTHIA 'W Monte Carlo generation (MC) 
followed by the simulation of the CDF II detector. QCD 
light jet production has a cross section about nine orders 
of magnitude greater than the theoretical expectation for 
the signal before the fo-tag requirement. Although this 
process generally does not have intrinsic ^t, mismea- 
sured jets can cause an imbalance in the total Et and 
fake the signal if one of the jets is misidentified as a 6 jet 
(mistag). Furthermore, 6 quark pair production yields 
h jets and I^t if one 6 quark undergoes a semi-leptonic 
decay. In both cases, the tends to be aligned parallel 
or anti-parallel to the first or second most energetic jet. 
We use MC simulation only to evaluate QCD processes 
which yield events with a 6 or c quark (heavy quark, h.q.) 
pair. This background is denoted as QCD h.q. in Table|l] 
The MC cross section for heavy quark jet production is 
scaled to the difference between the observed data and 
the predicted contribution of the other backgrounds. We 
determine the normalization factors, /cqcd, before the 
data are divided into control and signal regions to be 
1.30±0.04 and 1.47±0.07 for events with one and two h- 
tagged jets, respectively [T. The background due to light 
quarks being misidentified as b jets is denoted as mistags 



in Table [T] and is determined using a 6-jet misidentifica- 
tion rate that is measured from inclusive jet samples [8]. 
The ti cross section of 7.3 ± 0.8 pb is obtained from the 
combined CDF measurements [TD]. A small contribution 
by single top quark production is included, normalized 
with theoretical NLO production cross sections \n\ I12j . 
The electroweak backgrounds were generated at leading 
order with pythia and increased by a factor of 1.4 to 
account for higher order effects [T31 [H] . 

Simulated distributions of kinematic variables such as 
the jet Et, the dijet invariant mass, and the ifx have 
been found to be in agreement with observations in the 
control regions [7|. In Fig.[l] (a) and[l](b) we show the 
dijet invariant mass distributions when both jets are b- 
tagged in control regions CRl and CR2, respectively. 

We optimize the sensitivity to Higgs boson production 
by maximizing the ratio of iVniggs/V-^bck where A'niggs 
and -/Vbck are the number of predicted Higgs boson and 
background events in the signal region. The final se- 
lection requires A0(i?i^T,^T) > 0.8, ^t/Ht > 0.45, 
Et^i > 60 GeV , and > 70 GeV . The variable Ht 
{^t) is the scalar (the magnitude of the vectorial) sum of 
the two jet energies in the transverse plane. The number 
of background events after the final selection is shown in 
Table |l] For a Higgs boson of mass 115 GeV/c^ we ex- 
pect a total of 1.2 (0.4) events with one (two) 6-tagged 
jets [H]. 

The total systematic uncertainties on the background 
predictions are 17% and 19% in the single and double tag 
events, respectively [7|. One of the dominant correlated 
uncertainties for the background and signal predictions 
is due to the jet energy scale (JES) [6J. The JES un- 
certainty affects selection efficiencies non-uniformly over 
kinematic distributions leading to uncertainties in the 
distributions' shapes. The resulting systematic uncer- 
tainties are typically between 10% and 20% for multi-jet 
events depending on the kinematic selections, 15% for 
dibosons, 26% for W plus h.q., 17% for Z plus h.q., 1% 
for top quark, and 8% for signal events. The system- 
atic uncertainty on the ^-tagging efficiency is 4.3% for 
single and 8.6% for double 6-tagged events. There is a 
3% uncertainty on the trigger efficiency, 2% uncertainty 
on the signal and background acceptances due to uncer- 
tainties in the parton distribution functions, and 6% on 
luminosity. The dominant uncorrelated systematic un- 
certainties are due to the cross sections and the statistics 
of the simulated samples. We assign uncertainties of 11% 
to the top, 11.5% to the diboson, and 40% to the W plus 
h.q. and Z plus h.q. cross sections. Statistical uncer- 
tainties due to Monte Carlo sample sizes are 20%-32% 
on the multi-jet, 20% on the W plus h.q., and 11% on 
the Z plus h.q. predictions. Correlated and uncorrelated 
uncertainties are evaluated separately and combined in 
quadrature. 

The signal region is analyzed after all the background 



6 



TABLE I: Comparison of the total number of expected and observed single and double 6-tagged events in the control regions 
and in the optimized signal region. 
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See text 
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0.5 ±0.4 


18.3 ±8.1 
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FIG. 1: Dijet invariant mass distribution for double-tagged events in (a) CR 1 and (b) CR 2 



predictions and selection requirements were determined. 
The results are summarized in Table |l] Requiring one 
(two) 6-tagged jets we observe 268 (16) events, where 
248 ± 43 (14.4 ± 2.7) are expected. The dijet invariant 
mass distributions in the signal region when one and two 
jets are ^-tagged are shown in Fig. |2] 

Since no significant excess is observed, we compute 
95% C.L. upper limits for the Higgs boson production 
cross section times the branching fraction when the Higgs 
boson is produced in association with a, Z or W boson 
and decays to two b quarks while the Z and W yield large 
I^T and no identified leptons. The limits are computed 
using the Bayesian likelihood method pj] with flat prior 
probability for the signal cross section and Gaussian pri- 
ors for the uncertainties on acceptance and backgrounds. 
The likelihood quantifies the agreement between the dijet 
mass spectrum in the observed data and the sum of the 
modeled background processes and VH signal as a func- 
tion of the Higgs boson mass. We combine the channels, 
such as single 6-tag, double 6-tag, ZH, and WH by tak- 
ing the product of their likelihoods and simultaneously 
varying the correlated uncertainties. Table [H] summa- 



rizes the combined single and double 6-tag upper limits 
separately for ZH and WH , and the combined VH up- 
per limits. The result is also expressed as the ratio of 
the upper limit to the SM cross sections at NLO as a 
function of the Higgs boson mass. The observed limits 
agree with the expected ones at the 1 a level. 

In summary, we have performed a direct search for the 
SM Higgs boson decaying into b jet pairs using data with 
integrated luminosity of 0.97 ± 0.06 fb~^ accumulated in 
Run II of the CDF II detector. We observe no signif- 
icant excess over the background predicted by the SM 
and thus set 95% C.L. upper limits for the cross section 
of the Higgs boson production. The observed (expected) 
combined limit of the weak boson associated Higgs boson 
production cross sections set by this analysis for a Higgs 
boson mass of 115 GeV/c^ is 20.4 (14.2) times the SM 
prediction. These results significantly improve the com- 
bined limits on the allowed production rate for the SM 
Higgs boson at the Tevatron. 

We thank the Fermilab staff and the technical staffs 
of the participating institutions for their vital contribu- 
tions. This work was supported by the U.S. Department 
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FIG. 2: Dijet invariant mass in the signal region for (a) single-tagged and (b) double-tagged events. The dotted line overlaid 
on the background shows 10 times the expected yield of 115 GeV/c^ Higgs bosons. 



TABLE II: The expected cross section limits of the ZH and WH processes, and the expected and observed combined (VH) 
cross section limits when H — > bb. The last column gives the ratio of the observed limit with respect to the SM cross section. 
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